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ABSTRACT 

We have discovered an Ofpe/WN9 (WNll following Smith et al.) star in 
the Sculptor spiral galaxy NGC 300, the first object of this class found outside 
the Local Group, during a recent spectroscopic survey of blue supergiant stars 



^Based on observations obtained at the ESO Very Large Telescope [65.N-0389(B) and 67.D-0006(A)]. 
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obtained at the ESO VLT. The hght curve over a five-month period in late 
1999 displays a variability at the 0.1 mag level. The intermediate resolution 
spectra (3800-7200 A) show a very close resemblance to the Galactic LBV AG 
Car during minimum. We have performed a detailed non-LTE analysis of the 
stellar spectrum, and have derived a chemical abundance pattern which includes 
H, He, C, N, O, Al, Si and Fe, in addition to the stellar and wind parameters. The 
derived stellar properties and the He and N surface enrichments are consistent 
with those of other Local Group WNll stars in the literature, suggesting a similar 
quiescent or post-LBV evolutionary status. 

Subject headings: galaxies: individual (NGC 300) — galaxies: stellar content — 
stars: Wolf-Rayet 

1. Introduction 

With the new telescopes of the 8-10 meter class stellar astronomy is branching out 
beyond the Local Group. Ideal targets for our understanding of young stellar populations 
in distant galaxies are hot massive stars. These objects have strong stellar winds producing 
broad and easily detectable spectral features distributed over the whole wavelength range 
from the UV to the IR, and providing unique information on chemical composition, galactic 
evolution and extragalactic distances. With these perspectives in mind we have recently 
begun a systematic spectroscopic study of luminous blue stars in galaxies beyond the Local 
Group, and presented spectral classification and first quantitative results for A supergiants 
in NGC 3621 (6.7 Mpc, Bresohn et al. 2001) and NGC 300 (2.0 Mpc, Bresohn et al. 2002). 
Here we report on the discovery and detailed quantitative analysis of the first example of an 
Ofpe/WN9 star outside of the Local Group. We will present a detailed chemical abundance 
pattern - the first in a galaxy beyond the Local Group - together with stellar parameters 
and a determination of the stellar wind properties. 

The Ofpe/WN9 class was introduced to include objects which show in their spectra high 
excitation emission lines from Hell and Nlll, typical of Of stars, together with low excitation 
lines from He I and Nil, seen in late WN stars (Walborn 1982, Bohannan & Walborn 1989). 
Objects of this class have so far been identified in the Galaxy (possibly several stars in the 
Galactic center: Allen et al. 1990, Najarro et al. 1997, Figer et al. 1999), the LMC (ten 
stars: Bohannan & Walborn 1989), M33 (seven stars: Massey et al. 1996, Crowther et al. 
1997) and M31 (one star: Massey & Johnson 1998). The importance of Ofpe/WN9 stars as 
objects in a transitional stage of evolution between O and W-R stars has been recognized 
in the last decade, and a connection to the LBV class has been suggested, Ofpe/WN9 stars 
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being observed during a quiescent or post-LBV phase (Crowther & Smith 1997, Pasquah et 
al. 1997). Indeed, in at least a couple of instances Ofpe/WN9 stars have been observed to 
turn into LBVs (R127: Stahl et al. 1983; HDE 269582: Bohannan 1989). The LBV AG Car 
is also known to show an Ofpe/WN9-like spectrum during its hot phase at visual minimum 
(Stahl 1986). The discovery of ejected circumstellar nebulae, in some cases measured in a 
state of expansion, associated with some of the LMC Ofpe/WN9 stars by Nota et al. (1996) 
and Pasquali et al. (1999) brings forward strong evidence for the occurence of violent episodes 
of mass loss in these stars, similar to the shell-producing, eruptive outbursts of LBVs. 

Smith et al. (1994) revised and extended the classification of late WN stars to include 
lower excitation objects. In their scheme, stars showing spectra like AG Car at minimum, 
where no Nlll is detected, are reclassified as WNll. Given the spectral resemblance of the 
NGC 300 star we analyze here to AG Car at visual minimum (see § 4) we will adopt the 
WNll classification for the remainder of this Letter. We describe the observational data in 
§ 2, and the photometry in § 3. In § 4 we present our VLT spectra, together with the stellar 
and wind properties derived from a quantitative spectral analysis. 

2. Observations 

The data presented here are part of a spectroscopic survey of photometrically selected 
blue supergiants in the Sculptor spiral galaxy NGC 300 obtained at the VLT with the FORS 
multiobject spectrograph, and described in detail by Bresolin et al. (2002). This latter paper 
presents spectra obtained in September 2000 in the blue spectral region (~ 4000-5000 A) and 
used for the spectral classification of about 70 supergiants. The emission-line star analyzed 
here corresponds to star B-16 of the spectral catalog (see Table 2 and finding charts in 
the aforementioned paper), a rather isolated bright star apparently not associated with any 
prominent OB association, star cluster or nebulosity. The nearest OB association (from the 
catalog of Pietrzyhski et al. 2001) lies at a deprojected distance of ~450 pc. Given the 
estimated age of 3.8 Myr (§4.1) a speed on the order of 100 km/s would be required for the 
star to reach its current position from this association. This is much larger than the typical 
stellar random velocities, and comparable to that of Galactic runaway stars. To the best of 
our knowledge the star is not included in published listings of emission-line or W-R stars in 
NGC 300 (Schild & Testor 1992, Breysacher et al. 1997, and references therein), and we will 
refer to it as star B-16. 

In a more recent observing run (September 2001), in order to measure the mass-loss 
rates of blue supergiants in NGC 300 from the Ha line profile, we secured MOS spectra 
in the red with a 600 lines mm~^ grating and 1" slitlets, which, together with the older 
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spectra, provides us with complete coverage at i? ~ 1, 000 resolution of the 3800-7200 A 
wavelength range, except for a 50 A-wide gap centered at 5025 A. The total exposure time 
was 13,500 s for both the blue and the red spectra, while better seeing conditions favored 
the blue observations. The mass-loss rates we derive for the whole blue supergiant sample 
will be discussed in a forthcoming paper, and we concentrate here on the single B-16 star. 

3. Photometry 

As part of a multi-epoch photometric study of NGC 300 carried out for the discovery 
and monitoring of Cepheid variables (Pietrzyhski et al. 2001, see also Pietrzyhski et al. 
2002), BV photometry from ESO/MPl 2.2m telescope WFI images is available for B-16 
covering nearly a one-half year period. The resulting B and V light curves arc shown 
in Fig. 1. An apparently irregular variability, on the order of 0.1 mag (peak-to-peak), is 
detected in both bands, at a significant level above the typical photometric uncertainty of 
0.015 mag. The color index remains approximately constant around B — V — —0.07. While 
bona fide LBV variables show large amplitude variability on timescales of years (Stahl et 
al. 2001), partly characterizing the outburst phenomena in these objects, smaller- amplitude 
variability is observed in Ofpe/WN9 stars or LBVs during quiescence (e.g., Stahl et al. 1984, 
Sterken et al. 1997). Assuming the intrinsic color index from the stellar atmosphere models 
{Eb-v = 0.06), R = Ay / Eb^v = 3.1 and an average magnitude V = 19.00, we obtain from 
the adopted distance modulus for NGC 300, m — M = 26.53 (Freedman et al. 2001), an 
absolute magnitude My ^ —7.72 ± 0.13. 

4. Spectral analysis 

The blue and red portions of the rectified VLT optical spectrum are shown in Fig. 2. In 
the top panel we have also superimposed the spectrum of AG Car during minimum, taken 
from the Walborn & Fitzpatrick (2000) atlas, in order to illustrate the remarkable similarity 
between the two spectra. Line identification is provided for most of the recognizable features, 
which characterize B-16 as a WNll star. These include the hydrogen Balmer series fines 
(pure emission). He I lines (together with Hell A 4686) and mostly low-excitation metal fines 
such as Nil, Felll and Sill. P-Cygni profiles are seen in most of the He I fines, as well 
as in Si III AA 4552-4575, Nil AA 4601-4643 and Fe ill AA 5127,5156. The absence of N ill is a 
discriminant factor against an earlier (WN9-10) classification. The WNll classification is 
further confirmed by the observed equivalent widths of He I A 5876 (32 A) and He II A 4686 
(lA), compared to those given by Crowther & Smith (1997). 
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4.1. Atmosphere models 

For the quantitative analysis of the stellar spectrum we have used the iterative, non-LTE 
line blanketing method presented by Hillier & Miller (1998). The code solves the radiative 
transfer equation in the co-moving frame for the expanding atmospheres of early-type stars 
in spherical geometry, subject to the constraints of statistical and radiative equilibrium. 
Steady state is assumed, and the density structure is set by the mass-loss rate and the 
velocity field via the equation of continuity. The velocity field (HiUier 1989) is characterized 
by an isothermal effective scale height in the inner atmosphere, and becomes a (3 law in the 
wind (e.g., Lamers et al. 1996). Better fits are obtained if instead of the standard (3 law a 
combined 2-/3 law is used, where each value of (3 (given in Table 1) characterizes the shape 
of the velocity field in the inner and outer parts of the wind, respectively. We allow for 
the presence of clumping via a clumping law characterized by a volume filling factor /(r), 
so that the 'smooth' mass-loss rate, Ms, is related to the 'clumped' mass- loss rate, Mc, 
through Ms — Mcj j^^'^ ■ The model is then prescribed by the stellar radius i?*, the stellar 
luminosity L*, the mass-loss rate M, the velocity field i'(r), the volume fiUing factor /, and 
the abundances of the elements considered. The reader is referred to Hillier & Miller (1998, 
1999) for a detailed discussion of the code. 

The main stellar parameters of our best- fitting model are summarized in Table 1, 
whereas fits to some of the most important line diagnostics are illustrated in Fig. 3. As 
can be seen, the majority of the spectral features are well reproduced by the adopted model. 
The Fe III AA 5127,5156 lines, together with Hell A 4686, provide excellent constraints for the 
effective temperature. Using the strength of the electron scattering wings of the Balmer lines 
we derive a clumping factor / = 0.25. Thus the "smooth" mass loss rate would be twice 
the value displayed in Table 1. The parameters describing the stellar wind, i.e., the terminal 
velocity v^o — 325 kms~^, the mass-loss rate Mc — 4.6 x 10~^ Mq yr~^, and the wind perfor- 
mance number r\ — Mvoo/ (L^/c) = 0.84, lie in the range found for Local Group WNll stars 
(Crowther et al. 1997). We estimate the uncertainty in L and Mq/ f^^'^ to be approximately 
10% and 30%, respectively. As a further comparison. Fig. 4 shows the location of star B-16 
in the H-R diagram, together with additional WN9-11 stars in the Local Group (Crowther, 
Hillier & Smith 1995; Smith, Crowther & Wilhs 1995; Crowther et al. 1997), as well as LBVs 
at minimum in the Galaxy (P Cyg: Pauldrach & Puis 1990; AG Car: Smith et al. 1994) and 
in the LMC (R71: Lennon et al. 1994). 

Table 2 summarizes the model fractional mass abundances. Uncertainties in the latter 
range between 0.2 dex (He, N, Si and Fe) and 0.3 dex (C, O and Al). The chemistry of 
B-16 resembles that of other known Local Group WN9-11 stars. The severe depletion of 
hydrogen inferred from our analysis (mass fraction X = 0.27), and the correspondingly large 
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helium surface abundance (H/He = 1.5 by number), are fairly typical for late WN stars 
(H/He = 0.8-3 by number), and are indicative of heavy mass-loss stripping of the stellar 
outer layers during the post-main sequence evolution of massive stars. We find a total CNO 
abundance close to that of Galactic main sequence B stars (Gies & Lambert 1992, Kilian 
1992), consistent with a half solar chemical composition in these elements. This agrees with 
the empirical O abundance derived for an H ll region located 30" away (region 5 observed 
by Pagel et al. 1979). On the other hand, Fe has a roughly solar abundance, suggesting an 
a/Fe ratio different than in the Galaxy. In comparison with the Galactic B star abundance 
pattern, N is overabundant (by mass) by a factor of 10, while G and O are reduced by factors 
of ~2 and 6, respectively, although the latter abundances are highly uncertain. 

The current evolutionary status of B-16 can be estimated by means of the recent Geneva 
stellar tracks including rotation (Meynet & Maeder 2000) with an initial rotational velocity 
Vint — 300 kms~^ and Z = 0.02 (Fig. 4). The star's position in the H-R diagram approaches 
the 60 Mq track, and an initial mass of ~55 Mq can also be derived from the relation between 
stellar luminosity and initial mass for WNL stars given by Schacrcr & Maeder (1992). From 
the 60 Mq stellar model more closely matching the observed H and He surface abundances 
we infer an age of 3.8 Myr and a present mass of ~36 Mq. At this stage, the predicted N 
overabundance with respect to the ZAMS abundance is 12 (mostly occuring during the MS 
phase, as a consequence of rotational mixing), in good agreement with our finding. 

According to its He- and N-enriched surface chemistry, B-16 might be in a dormant or 
post-LBV phase of evolution. This is also indicated by the position in the H-R diagram, 
well above the Humphreys-Davidson limit (Humphreys & Davidson 1979), and intermediate 
between the LBVs AG Gar and P Gygni, which have comparable mass- loss rates to B-16, as 
well as similar wind performance numbers. The NGG 300 star shows higher He enrichment 
at the surface (n^^g/nij ^ 0.7, versus 0.4-0.5 for the LBVs), and a somewhat faster wind 
(325 versus 185-250 kms"^; Najarro, HiUier & Stahl 1997, Langer et al. 1994), so that a 
more advanced (post-LBV) state seems more likely. This picture is in agreement with the 
conclusions reached by Crowther, Hillier & Smith (1995) and Growther et al. (1997) regarding 
the close connection between WN9-11 stars and LBVs. 

An unexpected reward out of our spectroscopic survey of blue supergiants in NGG 300, 
the WNll star B-16 studied here is the first star in this galaxy (and beyond the Local 
Group) for which we have attempted a detailed quantitative analysis. In the near future the 
systematic study of the massive stellar population in NGG 300 and similar galaxies within 
a few Mpc, well within current observational capabilities, will certainly provide new insights 
into massive stellar evolution and stellar abundances. 
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Fig. 1. — B and V light curves for NGC 300 star B-16, spanning an approximately five 
month period between September 1999 and January 2000. The B — V color index is shown 
in the top panel, the V (full dots) and B (open triangles) magnitudes in the bottom panel. 
Photometric errors (~0.015 mag) are comparable to the size of the dots. 
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Fig. 2. — Rectified VLT/FORS optical spectra obtained at blue (top) and red (bottom) 
wavelengths. Line identification is provided in the top section of each panel. The comparison 
spectrum of AG Car at visual minimum (courtesy N. Walborn and E. Fitzpatrick) is shown 
as a dashed line superimposed on top of the blue spectrum of star B-16. 
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Fig. 3. — Different spectral diagnostics in the observed spectrum of B-16 compared to our 
best fit model (dashed line), convolved with a Gaussian profile to match the spectral resolu- 
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WNll stars (open circles, labeled) and WN9-10 stars (open triangles; from Crowther et 
al. 1995; Smith et al. 1995; Crowther et al. 1997), as well as LBVs at minimum (open stars, 
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line), and the theoretical Geneva stellar tracks including rotation {vini — 300 kms~^) at solar 
metallicity, for M = 40 Mq and 60 Mq (dotted lines). These tracks are plotted from the 
zero-age main sequence to the point where the fractional mass of hydrogen X ~ 0.20. 
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Table 1. Basic properties and model parameters 



R.A. (J2000) 00 54 44.64 

Dccl. (J2000) -37 42 39.13 

V (mean) 19.00 

B-V (mean) -0.07 

Mv -7.72 

Model fit parameters 

Eb-v 0.06 

\oz{U/LQr 5.94 

T, (K)'' 26000 

Teff (K)b 22000 

R^/Rq"^ 46 

R2/3/RQ' 64 

log(Mc) (M0 2/r-i) -4.34 

i^oo (fems-^) 325 

^1 6.0 

132 0.9 

/ 0.25 

McVoo/{L^/c) 0.84 

a.^ ~ 20 
V ~ 2/3 



Table 2. Model abundances 



Species 


Relative number 


Mass 


X/Xq 




fraction 


fraction 




H 


1.5E+00 


2.7E-01 


0.4 


He 


l.OE+00 


7.2E-01 


2.7 


C 


3.6E-04 


7.7E-04 


0.3 


N 


2.5E-03 


6.3E-03 


7.7 


O 


5.5E-04 


1.6E-03 


0.2 


Al 


5.0E-06 


2.4E-05 


0.4 


Si 


8.5E-05 


1.7E-04 


0.2 


Fe 


1.5E-04 


1.5E-03 


1.3 



